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Persistence in Random Walk in Composite Media 
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Abstract. - We consider a class of inhomogeneous media known as composite media that is 
often encountered in experimental sciences and investigate the persistence probability of a random 
walker in such a system. Analytical and numerical results for the crossover time scales has been 
obtained for a composite system with two homogeneous components and three homogeneous 
components respectively. 



The phenomenon of persistence in various stochastic 
processes has been well documented over the past decade 
[1]- [13]. Even the most simple of all stochastic processes, 
a random walker in a homogeneous and infinite media, ex- 
hibits the phenomenon of persistence and the non trivial 
persistence exponent 9 has the value 1/2 [9]. In an ex- 
perimental setup finite boundaries become important and 
we have recently investigated the effect of finite bound- 
aries on the survival probability of a random walker in an 
homogeneous system [16]. Its then a natural question to 
ask as to how the survival probability behaves for a ran- 
dom walker in an heterogeneous system. Although ran- 
dom walk in spatially disordered media has been already 
studied [17]- [19] and in few cases exact results are known 
for a similar quantity ,the first passage time [20]- [?], lit- 
tle is known about the persistence probability. We shall 
consider one such class of heterogeneous system which is 
known as composite media and is often encountered in ex- 
perimental science. A composite system essentially com- 
prises of segments of different homogeneous media which 
differ in their macroscopic properties, such as diffusion 
coefficients. Redner has already investigated the first pas- 
sage properties for a diffusive process in such a composite 
system [?] . He considers a linear chain of N blocks each of 
length Li having diffusivities The mean first passage 
time in such a system is calculated to be 



1 — 1 



L{Lj 



(1) 



While the first passage probability is simply the proba- 



bility that the particle escapes from one of the boundaries, 
the persistence probability is different and is defined as 
the probability that the random walker has not crossed 
the origin up to time t. If the system was homogeneous 
then the persistence probability of a random walker would 
be simply p(t) ~ t~ e with = 1/2. A composite system 
is slightly different in the sense that near the boundaries 
the difference in diffusivities tend to give a bias to random 
walker. 

The simplest of composite media that can be con- 
structed is the one with two homogeneous segments with 
different diffusivities. We shall first derive the result for 
such a system and later generalize this result for different 
types of composite media. 

Consider two homogeneous media of diffusivities D\, 
henceforth called medium 1 and D2, medium 2. A slab 
of medium 1 is placed between —L and +L and the rest 
of the space is filled with medium 2 as shown in Fig. 1. 
Since the diffusion coefficients are different, it follows that 
the stochastic noise correlator will be different and in par- 
ticular they are 



(2) 



where i,j are the medium indices 1 and 2 and t\ > ti. 
For a random walker the the probability that the walker 
is at (x,t) starting from (xo,0) simply obeys the diffusion 
equation in two segments as 



W email: tpdc2@mahendra.iacs. res. in 



dP _ d 2 P 

dt 1 dx 2 

dP _ d 2 P 

~dt ~ 2 ~dx Y 



L<x>L 
L > x > L 



(3a) 
(3b) 
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Fig. 1: Arrangement of two homogeneous media. 

The exact dynamics of the problem can solved by con- 
sidering the Laplace transform of Eq.(3a) and Eq.(3b) in 
which case the solution to the equations become 



(4a) 



P(x,s) = Aicxp ( -J — \x\ 



L < x > L 



+^ 2 cxp ( J — \x\ 



P(x,s) = A 3 cxp \-\I-qA x 



(4b) 



The coefficients A\,A 2 ,A% arc found from the bound- 
ary conditions that the probability P(x,t) and the current 

dx 



-D^f- is continuous across the boundary. Finally, the 



third unknown coefficient is found from the normalization 
of the probability. The resulting expression, however, is 
complicated and it is difficult to extract any information 
from it. 

We, instead, take a different approach to derive our re- 
sult. The equation of motion for the random walker is not 
changed in spite of the heterogeneity of the system and is 
simply 



da; 

dt 
dx 

lit 



for 
for 



L<x<L, 
- L> x > L. 



(5a) 
(5b) 



Let the time required for a random walker to reach the 
boundary x = ±L be r. In which case we can write down 
the solution for the equation of motion in the two different 
regions. For — L < x < L the solution is simply 



x(t) = f dt' m (t'), 
Jo 



(6) 



whereas for — L > x > L, when the particle is in medium 
2 the solution for x{t) becomes 



x(t) 



£ dt' Vl (t') + J* dt' m (t'). (7) 



Eq.(7) simply states the fact the random walker has spent 
time t in medium 1 and the rest of the time in medium 2. 
Both Eq.(7) and Eq.(8) are valid when the walker is deep 
inside medium 1 or medium 2, since none of the equations 
considers the hopping across the boundary. When deep 
inside either media the multiple boundary crossings are 



rare events whereas when the walker is near the boundary 
multiple crossings arc frequent and it is due to these mul- 
tiple crossing events the crossover is not sharp and there 
will be two crossover time scales in the problem. 

The correlation (x(ti)x(t 2 )} can now be worked out 
carefully. First consider the case — L < x(ti) < L, 
—L < x(t 2 ) < L and t\ > t 2 . In this case the correla- 
tor becomes 

(x{h)x{t 2 )) = / 1 dt[ f 2 dt' 2 (m(t[)v(t 2 )) = 2D Y t 2 

Jo Jo 

(8) 

where we have used Eq.(2) for the noise correlator. If, 
however, — L > x(ti) > L, —L < x(t 2 ) < L and t\ > t 2 
then we have 



{x{h)x{t 2 )) = ( 



/ dt 2?7l (t 2 ) 
Jo 



) 



(9) 



Since {rii{t)rj 2 {t)) = 0, the above expression simplifies to 

(x(t 1 )x(t 2 ))= f dt[ [ 2 dt' 2 ( m (t' 2 ) m (t[)) (10) 
Jo Jo 

As r > t 2 the correlation {x{t\)x(t 2 )) becomes 



(z(fi)a:(t 2 )> = 2£) 1 i 2 



(11) 



Finally, for — L > x(ti) > L and — L > x(t 2 ) > L,ti> t 2 , 
the correlator becomes 



{x{h)x{t 2 )) = ( 



J dt' 1 r ?1 (i' 1 ) + ^ dtiffe(ti) 
J dt' 2 r h {t' 2 ) + J dt' 2V2 {t' 2 ) 



> (12) 



Since the cross-correlation of the noise is zero we arrive at 

{x(ti)x(t 2 )) = f dti / dt' 2 ( m (t[) m (t' 2 )) 
Jo Jo 

+ j* dt[ y* dt'Mt'MQ) (13) 

The first term in Eq.(13) is easy to calculate and gives us 

/ dti/ dt' 2 ( m (t' 1 ) m (t' 2 ) = 2D 1 T (14) 
Jo Jo 



To evaluate the second term we make a transformation of 
variable t" = t' — r and we have 



pt\—T rt^—T 
/ < / , 

Jo Jo 



dt'l2D 2 5{t'l - 4) = 2D 2 (t 2 - t) (15) 



Hence, the complete correlator is 

(x{h)x{t 2 )) = 2(Di - D 2 )t + 2D 2 t 2 (16) 
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Fig. 2: Plot of {x 2 (t)} vs time in log-log scale for L = 20, 2Di = 
20, 2Z?2 = 2. The crossover time scales and the crossover 
regime is indicated in the figure. The square points are actual 
data from numerical simulation and the circular points are fit 
of Eq.(17) 



Fig. 3: Plot of (x 2 (t)) vs time in log-log scale for L — 50, 2Di — 
50, 2D2 = 5. The crossover time scales and the crossover 
regime is indicated in the figure. The square points are actual 
data from numerical simulation and the circular points are fit 
of Eq.(17) 



Of all the quantities in Eq.(16) the only unknown is r. is then given by 
Since Eq.(16) gives us noise averaged quantities we might 
as well replace r by an average value, which is simply 
r = 275-, the average time for a random walker to reach 
x = ±L. 

It is clear from Eq.(ll) and Eq.(16) that there are 
two relevant time scales in the problem. The first one 

is Ti = ( Di d^ 2 ) 27J7 wnereas the second one is r 2 = 



{X(h)X(t 2 )) = for t 2 <n 

ii 



l[3L 2 + 2D 2 t 2 
[3L 2 + 2D 2 h 



for t 2 > t 2 , 



(18) 



It is between these two time scales when w ^h P 



\ £>i J 2D 2 

the random walker undertakes multiple hoppings across 
the boundary, as a result of which the temporal regime 
T\ < t < t 2 gives the crossover region in the system. 
The complete correlator now takes the form 



2D x t 2 



for t 2 < Ti, 



Di-Ds 
Di 



L 2 + 2D 2 t 2 for t 2 >r 2l 
(17) 



In the limit D x = D 2 = D we get the correct result for 
a homogeneous medium. 

A numerical simulation of the system for two different 
values of L and two different sets of D\ and D 2 has been 
done. Simulation result for the mean square displacement 
(x 2 (t)) is shown in Fig. 2 and Fig. 3. Configuration averag- 
ing of 10 4 has been done for both the systems to obtain 
the result. 

To calculate the survival probability in the two regimes 
we use Eq.(17) and with suitable transformations both in 
x and t convert the process to a Gaussian stationary pro- 
cess. Define a normalized variable X(t) = x(t)/ -J (x 2 (t)). 
The correlator in this normalized variable, (X(ti)X(t 2 )), 



D1-D2 



. For t < Ti we define the usual trans- 



formation in time, T = ln< and the two time correlation 
function in the new time variable becomes 



(X(T0X(T 2 )) = e 



-l/2(7i-T 2 ) 



(19) 



and the survival probability p(t) for this temporal regime, 
in real time, is then 

p(t)^t-^ 2 . (20) 
For t > t 2 we define a new time variable T as 



(3L 2 + 2D 2 t. 



(21) 



The correlation function (X(T\)X(T 2 )) takes the form of 
Eq(19), except that the time transformations are differ- 
ent. Since the process is a Gaussian stationary process 
and the correlator is exponentially decaying, the survival 
probability in the new time variable, P(T), is 



P{T) = e 



= p -T/2 



(22) 



In real time the survival probability p(t) takes the form 

1 



P(*) 



y//3L 2 + 2D 2 t' 



(23) 
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Fig. 4: Plot of survival probability p(t) vs time in log-log scale 
for L — 20, 2D 1 — 20, 2L>2 = 2. The crossover time scales and 
the crossover regime is indicated in the figure. The circular 
points are actual data from numerical simulation and the solid 
lines are fit of Eq.(21) and Eq.(??) 



Fig. 5: Plot of survival probability p(t) vs time in log-log scale 
for L — 30, 2Di — 30, 2Z?2 = 3. The crossover time scales and 
the crossover regime is indicated in the figure. The circular 
points are actual data from numerical simulation and the solid 
lines are fit of Eq.(21) and Eq.(??) 



A plot of the survival probability, Eq.(20) and Eq.(23) 
for the two time regimes is shown in Fig. 4 and Fig. 5. 
The crossover timescales T\, t 2 and the crossover regimes 
are also indicated in the figures. Configuration averaging 
of 10 6 has been done to obtain the numerical results of 
Fig. 4 and Fig 5. Theoretical and numerical values of t\ 
and t 2 are presented in Table A for two different set of 
values of L, D\ and D 2 . 



Table A 



Parameter Values 


T th 




~t/i 


nu 
T 2 


L = 20 
2L>i = 20 
2D 2 = 2.0 


18 


17.161 


180 


179.988 


L = 30 
2L>i = 30 
2D 2 = 3.0 


27 


27.228 


270 


271.899 



Survival probability for three media. 

In this section we consider a composite system that is 
made of three homogeneous media with diffusivities D\, 
D 2 and D 3 . The medium with diffusivity Di is placed 
between ±Li while the second medium with diffusivity D 2 
is placed symmetrically between — (L 2 + L\) < x < —L\ 
and Li < x < (Li + L 2 ) as shown in the figure. For a 
random walker in region I the average time to reach the 

L 2 

boundary ±Li is t = jjj^ . When the random walker is 
in region II the average time to cross a region the length 
L 2 is once again r' = Thus, for t < r the particle 

spends its time in region I, for r < t < t' the walker is in 
region II while for t > t' the walker escapes to region III. 



The equation of motion in all the three region are 
dx 

~r = Vi(t) Wltn i=l,2,3 for three media, (24) 
at 

with the noise correlator 

{r li (t)i lj {t'))=2D i S ij 8{t-t l ), (25) 

where i, j are the medium indices running from 1 to 3. 

The solutions to Eq(24) for the three regions are respec- 
tively 

x(t) = [ dt'mit') for -Li < x{t) < L x (26) 
Jo 

x(t) = J dt' m (t') + J dt' m (t') (27) 
for Li < x(t) < L 2 and — L 2 < x(t) < —L\, 

x(t)= f dt' m (t')+ f + dt' m {t')+ f dt' m {t') 

JO Jt Jt+t' 

for x{t) > L 2 and x{t) < —L 2 . 

(28) 

The two time correlation function (x(ti)x(t 2 )} can be 
worked out carefully and for both x{ti) and x(t 2 ) lying in 
region I, with t\ > t 2 is 

(xihMh)) = 2D x t 2 (29) 

For x(ti) and x{t 2 ) lying in region II, {x{t\)x{t 2 )) takes 
the form 

(x(h)x{t 2 )) =l3L\ + 2D 2 t 2 (30) 
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Fig. 6: Arrangement of three homogeneous media. 



while for x(t\) and x(t\) lying in region III, using the 
fact that the cross correlations of the noise is zero, the 
correlator becomes 



(x(h)x(t 2 )} 



+ 



J T-\-t' J T 





r+r 



d*a<»7i(*i)»7i(*2)) 
dt'Mt'Mt'2)) 

dt' 2 {m(t'Mt' 2 ))- 



(31) 



The first integral is simply 2D\t. The second integral is 
performed by making use of the transformation t" =t' — r 
and the integral reduces to 



f f 

Jo Jo 



dt' 2 (V2(t'Mt' 2 )) = 2D 2 r' 



while for the third integral we use the transformation 
t" = t' — (r + t') and the integral is evaluated to be 
2D 3 (t 2 — t — t'). Hence the correlator becomes 



(x{t{]x(t 2 )) 



2D iT + 2D 2 t' + 2D 3 (t 2 
{3 x L\+{3 2 L\ + 2D 3 t 2 , 



r') 



(32) 



with ft = (Di - D 3 )/D 1 and (3 2 = (D 2 - D 3 )/D 2 . It is 
clear from Eq.(29), Eq.(30) and Eq.(32) that there four 
relevant time scales in the problem. The first one is obvi- 

L 2 L 2 

ously n = /?2d7- The second one is r 2 = P 2 jj^- The 
temporal regime n < t < t 2 represents the crossover 
regime from region I to region II, when the walker feels 
the effect of the inhomogeneity. Similarly, the third time 
scale is t 3 = i^j^{(3 1 L\+ (3 2 L\) and the fourth time scale is 
T4 = j^(PiLl+f3 2 L 2 ). T3 < t < T4 is the crossover regime 
from region II to region III and it is during this time when 
the walker spends most of its time near the boundary of 
region II and region III. Thus the proper time scales for 
which Eq.(29), Eq.(30) and Eq.(32) are valid are respec- 
tively < t < n, t 2 < t < T3 and t > T4 while the time 
intervals n < t < t 2 and T3 < t < T4 represents the two 
crossover regimes. 



The mean square displacement (x 2 (t)} is then 

(x 2 (t)) = 2D x t forO<i<Ti 

= @L\ + 2D 2 t 2 for t 2 <Kt 3 

= p 1 L\ + l3 2 L\ + 2B ? ,t 
for t > r 4 . (33) 

Note that for D 2 = D 3 we recover the first case, that 
is a composite media with two homogeneous components 
while for D\ = D 2 = D3 = D we recover the case for a 
homogeneous system. 

To obtain the survival probability from Eq.(29), Eq.(30) 
and Eq.(32) we follow the usual procedure of defining suit- 
able transformations in space and time as in the earlier 
section. Thus, we define a normalized variable X(t) as 
X(t) — , x and the correlator in the normalized vari- 

able becomes 



<*(*l)*(*2)> 



'j- for < t < ti 



lpL\ + 2D 2 t 2 
0L{ + 2D 2 h 



for t 2 < t < r 3 



\(j 1 L\^(i 2 L\ + 2B ?> t 2 
fcLl + d 2 L\ + 2D 3 h 



for t > r 4 
(34) 



The time transformation t — ► T for the three different 
regimes are defined in the following way 

e T = t for < t < t\ 

= f3L\ + 2D 2 t for r 2 < t < r 3 
= f3iL\ + j3 2 Lj + 2D 3 t for t > r 4 



and the correlator (X(T\)X(T 2 )) becomes 



(X(T!)X(T 2 ) 



e 2 



|(Ti-T 2 ) 



(35) 



(36) 



for all the three temporal regimes, the difference being that 
the time transformations arc different in the three regimes. 
The process is now a Gaussian stationary process. 

Since the correlator is exponentially decaying, the sur- 
vival probability in the transformed time variable is simply 



P(T) 



-T/2 



(37) 



In real time, using Eq.(37), the survival probability be- 
comes 



P(t) 



r 1/2 for < t < ti 
1 



^{3L\ + 2D 2 t 
1 



for t 2 < t < t 3 



(38) 



for t > T4. 
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Fig. 7: Plot of survival probability p(t) vs time in log-log scale 
for Li = 60, L 2 = 360, 2£>i = 200, 2D 2 = 20 and 2D A = 2. 
The crossover time scales and the crossover regime is indicated 
in the figure. The circular points are actual data from numer- 
ical simulation and the solid lines are fit of Eq.(??). 



Fig. 8: Plot of survival probability p(t) vs time in log-log scale 
for Li = 30, L 2 = 210, 2L>i = 100, 2D 2 = 10 and 2D Z = 5. 
The crossover time scales and the crossover regime is indicated 
in the figure. The circular points are actual data from numer- 
ical simulation and the solid lines are fit of Eq.(??). 



A plot of the survival probability for two different set of 
values of L\, L2, D\, D2 and D3 is shown in Fig. 5 and Fig. 
6. Configuration averaging of 10 6 has been done to obtain 
the numerical results of Fig. 5 and Fig. 6. Theoretical 
and numerical values of the time scales are presented in 
Table D. 
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Table D 



Parameter Values 


Time Scales 


Li = 30 
L 2 = 210 
2D 1 = 100 
2D 2 = 10, 2D 3 = 5 


t\ h = 8.1, t[ iu = 8.68 


r\ h = 81, r 2 " u = 80.4393 


rf = 1705.5, r 3 ntl = 1702.58 


T\ h = 3411, tI u = 3421.27 


Li = 60 
L 2 = 360 
2L>i = 200 
2D 2 = 20, 2D 3 = 2 


r{ h = 16.2, rf" = 16.01 


r\ h = 162, t 2 " u = 162.096 


rf = 4228.2, t 3 " u = 4245 


r\ h = 42282, r] m = 42329.8 



To conclude, we have investigated the phenomenon of 
persistence for the case of a random walker in a compos- 
ite media with two and three homogeneous components. 
We have presented a very simplified theory to explain the 
survival probability of a random walker in such inhomoge- 
neous systems. For the two component system, analytical 
results show that there are two relevant time scales in 
the problem and this time interval is the crossover regime 
for the problem. Similarly, for the three component sys- 
tems there are four relevant time scales and two crossover 
regimes. The fact that the crossover regimes are not sharp 
is due to the multiple hoppings that a random walker un- 
dergoes near the boundary. 
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